Non-universality of Anderson localization in short-range correlated disorder 
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We provide an analytic theory of Anderson localization on a lattice with a weak short-range correlated disor- 
dered potential. Contrary to the general belief we demonstrate that already next-neighbor statistical correlations 
in the potential can give rise to strong anomalies in the localization length and the density of states, and to the 
complete violation of single parameter scaling. Such anomalies originate in additional symmetries of the lattice 
model in the limit of weak disorder. The results of numerical simulations are in full agreement with our theory, 
with no adjustable parameters. 
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It is customary to assume that a wave completely looses 
its phase memory when reflected from a weak disordered po- 
tential flllill. Many powerful analytical frameworks, 
such as the single-parameter scaling theory |5], the DMPK 
approach |6, 7], or ID non-linear er-model |8] are based, im- 
plicitly or explicitly, on the reflection phase randomization 
(RPR). The RPR stands behind the notion of the standard 
universality classes in the random matrix description of non- 
interacting disordered wires |9] and leads to the independence 
of the mean density of states on disorder strength. 

The violation of the RPR is responsible for many non- 
universal effects. One example is the tight-binding model with 
hopping disorder, where the density of states diverges at the 
band center for an odd number of coupled chains 1 10], while 
it vanishes for an even number 0. The breaking of single 
parameter scaling and RPR by hopping disorder and other de- 
viations from universality were studied recently in great detail 
00- 

A partial violation of the RPR can be induced by a disor- 
dered on-site potential alone. This happens when the potential 
can no longer be regarded as weak, e.g., at the edges of elec- 
tronic conduction (or photonic transmission) bands 0[3> 
or in the situation when the probability density of the potential 
has power law tails 1 16]. In the case of a one-dimensional lat- 
tice with white-noise potential the RPR is partially violated at 
the band center, leading to the Kappus-Wegner anomaly char- 
acterized by a 9 % increase of the localization length ll7U18ll . 

In this paper we demonstrate that the RPR can be broken 
far more dramatically if the disordered potential is short-range 
correlated. (By the short range we mean a finite correlation ra- 
dius which is much smaller than the localization length.) The 
lack of RPR is accompanied by anomalies in the localization 
length (which can sharply increase or decrease), and in other 
quantities such as the delay time or the density of states. In 
contrast to the general belief even next-neighbor sta- 

tistical correlations in the potential can lead to a severe vio- 
lation of the single parameter scaling. In brief we distinguish 
two different effects of the correlations: (i) the system retains 
its universal properties with a renormalized localization length 
12111 : (ii) the universality is broken in the vicinity of specific 
spectral points; the RPR and the single parameter scaling are 



FIG. 1: A lack of phase randomization leads to a distinctive band 
center anomaly in the conductance g n through a disordered wire Q 
of length n — 12000 with next-neighbor correlations I17i of the po- 
tential and a = 1/150 (RPR localization length ^ = 2/ a = 300). 
The main panel shows the mean logarithm ci = — (l/2n)(ln g„) 
and its variance C2 = (l/4n) Varlng n from numerical simulations 
(data points) and the analytical expressions I21l24ai (solid lines), 
scaled to The inset shows the distribution of the reflection phase 
P(a) at the band center, calculated numerically as well as analyti- 
cally (HJ. 



violated; the density of states shows an anomaly, which de- 
pends on disorder strength. 

We consider the one-dimensional Anderson model 



(1) 



with a weak disordered potential U n , (U n ) = with \U n \ <C 
1. 

The correlation function (U n U m ) — a(n,m) is assumed 
to be invariant under a finite translational shift a(n + Q, m + 
Q) = <r(n,m). The minimal period Q does not need to be 
identical with the lattice constant. This accommodates the 
cases of structural and chemical disorder, or carefully engi- 
neered disorder, e.g., with masks. For Q = 1 the correlation 
function is generally written as 



(U n U n 



(2) 
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The localization length £ is accessible via the dimension- 
less conductance (transmission probability) g n of the system 
of length n, = — lim n _ >00 (l/2ri)(liig , n ). In order to 
quantify deviations from universality we consider the com- 
plete statistics of the conductance fluctuations obtained from 
the generating function 
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X (/i)= lim - ln(g-^ 2 ) = 



fi > 0. (3) 



to the second order in the potential. Similarly, the increment 

SX = X n+q - A„ is given by 
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SX n =J2<+s(a + 2ks) 



(9) 
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-<$> n+s {a + 2fcsX +m (a + 2km). 



The localization length is given by £ = q 1 . In single- 
parameter scaling |5], = c\, and c s — for s > 3, corre- 
sponding to a log-normal distribution of g n . In this paper we 
concentrate on the first two coefficients c i and c? . We base 
our analysis on the exact phase formalism I22L 12311 , which we 
extend to the case of correlated disorder. 

For the sake of definiteness assume the wave is reflected 
from the right end of a system of length n ^> £, with reflec- 
tion amplitude r n = — g n exp(ia„). The reflection phase 
a n 6 (0, 2tt) is related to the wave function by 
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E 



while the conductance g n is obtained from 



-2 cos k, 



(4) 



£ |. (5) 



Up to the second order in disorder strength, the statistical evo- 
lution of the phase and conductance with increasing system 
size is described by the recursion relations 

a n+ i-a n = 2k - 2K n (a n ), (6a) 
A„+i-A„ = K' n (a n ) + (K' n (a n )) 2 , (6b) 

where the prime stands for the derivative with respect to a n 
and the function K n (a) is given by 

K n {a) = - ^] , , $„(«) = — (1 + cosa), (7) 

with the group velocity v = \dE/dk\. 

There exsist two major obstacles in the derivation of 
the Fokker-Planck equation for the joint probability density 
P n (a,X) from Eqs. <6al6bt . First, the wave-number k on 
the right-hand side of Eq. do*at is not small. Second, the 
values of K n are correlated at different sites. These diffi- 
culties can be circumvented by monitoring the variables a n 
and A n with a step of q sites ('24]. We parameterize the en- 
ergy E = — 2 cos np/q + e with small e -C 1 and integer 
p, and choose q to be much larger than the correlation radius 
of the potential. Thus, the change of the phase over q sites 
5a = a n+q — a n is small, 
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-<f>„ +s (a + 2fcs)<J>^ +m (a + 2km), 



In the limit o weak disorder the recurrent relations Eqs. ( I9I8> 
lead to the Fokker-Planck equation for the joint probability 
density P n (a, X) 



9P n d , s 1 d 2 , s 
~dn = ~~da ^ + 2~dcp ^ 



(10) 



-^■Fi(a)P n + \^D x {a)P n + ~D ltl in )/>„. 
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The drift and diffusion coefficients, which determine the 
phase distribution P n (a), are 

1 1 °° 

F(a) = -(San), D(a) = - (5a n Sa n+mq ). (11) 

m— — oo 

In the case of the correlated disoder (0 with Q = 1, these 
coefficients are related by F = 2e/w + (l/A)dD/da. The 
other coefficients in Eq. fllOI > are given by 

1 1 °° 

Pi (a) = -(<5A„), Pi (a) = - £ (<5A„<5A„ +m9 ). 
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An (a) 



(<5a n <5A n+m< j) 



(12) 



Thus, we have derived Fokker-Planck equations i ll Oi l, which 
describe the system in the vicinity of a given rational energy 
E = —2cosTrp/q with q > 2. The potential fluctuations 
are assumed to be restricted within the conduction band, so 
that \E ± U n \ < 2. In particular Eq. dlOi does not apply in 
the vicinity of the band edge, where any fluctuation is strong. 
The latter case has to be treated separately. The effect of di- 
chotomic correlated disorder near the band edge was studied 
inRef. El. 

The equation for the phase distribution function P n (a) is 
readily obtained by integrating Eq. dlOt over the variable 
A. There exists, however, no general way to derive a sim- 
ilar equation for the distribution P„(A). The calculation of 
the generating function x(m) hence requires the analysis of 
the full density P„(a, A). Such analysis is greatly simpli- 
fied for RPR, which implies the factorization P n (a, X) — 
(2tt)^ 1 P„(A) . In this case a closed equation for P n (A) can 
be derived by averaging Eq. JTol over the phase 



dPn(X) 

dn 



i dP n (x) , i d 2 p n {x) 



^ ox 2^ dx 2 



(13a) 



da 
2^ 



Pi (a) 



da 
2^ 



Pi(a), (13b) 



3 



where £^ is the localization length in the presence of RPR. 
The last equality in Eq. dl3b> follows directly from Eq. ( I6bt . 
which implies that the first and second moment of the incre- 
ment of A are equal if the phase is randomized. The generating 
function calculated from Eq. (II 3bi has the parabolic shape 
X(m) — Hence, RPR implies single-parameter 

scaling in the localized regime, even in the presence of corre- 
lated disorder [this statement holds also for hopping disorder, 
which modifies the expression for K n , while Eqs. ( I6al6bi re- 
tain their form.] 

It is instructive to calculate £^ for correlations of the type 
(0. Taking the integrals in Eq. (II 3bi we reproduce the result 
ofRef. Hi 
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2v 2 ^ 



(14) 



It is important to remember, however, that ^ = £ only if RPR 
holds. This brings us to the question: What are the conditions 
for RPR, and what are the implications for the localized wave 
functions when RPR does not occur? 

As a first example to illustrate the effect of short-range cor- 
relations on the reflection phase statistics, we consider dis- 
order of the type 10 with next-neighbor correlations only, 
Cs>3 = 0. We let q = 2p, q ^S> 1 and obtain from Eq. 



D(a) = 2<jq — (ctq — 2<7i) sin 2 a. 



(15) 



The Fokker-Planck equation for the stationary phase distribu- 
tion P(a) is simplified to 



6 da 



ld_ 

2da 



d 

d—Vdp = 0. 

oa 



(16) 



Its solution at the band center e = E = is given by P oc 
[£>(a)] _1 / 2 . Note, that the solution slightly deviates from the 
constant even for the white-noise potential (cii = 0), which is 
the source of the so-called Kappus-Wegner anomaly 1H4I ITtI 
Hill . The next-neighbor statistical correlations can induce far 
stronger deviations from RPR. Indeed, for a = ay = — <tq/2, 



(U n U n >) = 2a8 nn > - crS mv 



±1. 



(17) 



the solution to Eq. (1161 is singular at e = 0, because the func- 
tion D(a) develops a zero for a — ±ir/2. 

At the level of Eq. dl6l the situation is analogous to the 
band center Dyson singularity [ 10] in the presence of hopping 
disorder. Hence, the correlations turned a weak anomaly into 
a strong anomaly. While the usual Dyson singularity appears 
as a consequence of an exact "crural" symmetry of the wire, 
such symmetry is limited to the first two orders in disorder 
strength for the correlated disorder dl7> . Using the decompo- 
sition U n = u n — u„_i with (u n u n i) — (j5 nn i, and taking 
also the fourth order terms (u^) = rj 2 in the potential into 
account, we derive at the band center the regularized phase 
distribution 



P E =o(u) oc (cos 2 a + (if - a 2 )/Aa) 



-1/2 





(E-\)/2a 

FIG. 2: The mean logarithm of the conductance and its variance are 
calculated numerically for the system Q of the large length n. The 
coefficients ci = — (l/2n)(ln<jr n ) and C2 = (l/4n) Varlng™ are 
plotted as the function of energy near the quarter band E = 1. The 
disorder is generated by U n — 2u„ — u„+i — u n -i (top panel) 
and U n = 2u n +i — u n — u n -i (bottom panel), where the random 
numbers u n = unless n is a multiple of 3 and {u3 m u im i) — 
a8 mm i, a — 1/150. The insets show the distribution function of 
the reflection phase. The solid lines are obtained from the solution 
of Eq. i 101 with p = — 1, q = 3. The RPR localization length ^ is 
obtained from Eq. I13bi . 



which becomes more singular for lower disorder strength. The 
probability density Pe=o(c() acquires a non-universal depen- 
dence on the shape of the distribution function of the poten- 
tial via the relation between its fourth and second moment. 
Slightly away from the band center, e ^> \J oirj 1 — a 2 ), the 
fourth-order terms play no role and the phase distribution is 
described by the solution of Eq. H6\ with D = 4cr cos 2 a, 



poo e -ey/2a 

P(a) oc / dy ■. 

Jo v y 2 cos 2 a + y sin 2a + 1 



(19) 



(18) 



The lack of RPR is necessarily reflected in an anomaly 
of the density of states, due to the node-counting theorem 
ll26H27ll : this also implies a different statistics of time delays 
t = da/dE, for which the increments are obtained by differ- 
entiating Eq. l|8) with respect to e. 

We now explore the implications for the transmission prop- 
erties of the system. For the specific correlations ( 1171 . the 



4 



Fokker-Planck equation dlOi simplifies near E = by 



dP n 

dn 



dP n 
da 



2a 



d 1 . 8 

■ 7 — cos a + — sma— 
da 2 oX 



P n . (20) 



The coefficient c\ is obtained directly by averaging the ex- 
pression with the stationary phase distribution d!9l >. The 
result is 



I-2TT 

ci = a da P(a) cos 2 a 
Jo 



e T K 1 (ie/2a) 
- — 1m - ; — - . 

2 K {ie/2a)' 



(21) 



where K u (x) is the modified Bessel function. The singularity 
in Eq. J2 It in the limit e — > is again regularized in the fourth 
order of the potential. In the weak disorder limit we find at the 
band center £/£;0 = In8y/cr/(r) 2 — a 2 ). 

From Eq. (I20> we can determine all coefficients c s recur- 
sively double Laplace transform 



In conclusion we show that the ID Anderson model with a 
weak short-range correlated potential may demonstrate strong 
anomalies near specific spectral points. Such anomalies can 
be used to strongly modify the transmission properties of elec- 
tronic wires and photonic wave guides in very small energy 
windows provided the phase coherence preserved over a long 
distance. In these windows the reflected wave develops a pre- 
ferred phase relation with the incoming wave. These proper- 
ties indicate that disorder correlations may be favorably em- 
ployed in the design of photonic or electronic filters. 

We gratefully acknowledge discussions with Piet Brouwer 
and the warm hospitality in the Lorentz Institute in Leiden, 
where part of this work was done. 



P x (a,fx) 



E 



dX^~ x ^ )n P n (a,X), (22) 



reduces Eq. (I20> to the eigenvalue problem 

-I 2 



P_Px 
da 



2a 



d a . 

—- cos a sin a 

da 2 



P X =X(M)P X - ( 23 ) 



The generating function xil 1 ) can be obtained perturbatively 
in /j, taking the solution P(a) of Eq. dl9l as zero approxima- 
tion. In particular, the variance c 2 is given by 



(-2 

G(a) 
r 



a — ci — / da (ci — cr cos 2 a)G(a), (24a) 
Jo 



d_ 

da 
d 



C^ 1 (ci -a cos 2 a+a-^- sin 2a)P(a), (24b) 
e d d 



— — + — cosa— cosa. 
2a da da da 



In order to illustrate our predictions, we compare them in 
Fig-Eto the results of numerical simulations. The numerical 
results agree with the theory, without any adjustable parame- 
ter. The main panel shows the energy dependence of c\ and 
C2 near the band center [Eqs. (!21l24aH . which clearly deviates 
from the single-parameter scaling prediction ci = c 2 = £7 . 
The inset shows the phase distribution dl 81 at the band center. 

In general, RPR is completely violated if the diffusion co- 
efficient D(a) has zeros as the function of a. The reflected 
wave then mantains a strict phase relation with the incoming 
wave. For the correlated disorder of the type (|2) such relation 
(phase selection) can only occur at the band center. The re- 
striction is lifted for correlations with Q ^ 1. In Fig. (|2j we 
provide examples of a quarter band anomaly, caused by a dis- 
order correlations with Q = 3. Following this recipe, strong 
anomalies can be produced in a vicinity of arbitrary rational 
values of the wave length 2q/p by a suitably correlated weak 
disorder potential with Q = q. This is in striking contrast to 
the case of white-noise disorder, which produces only a weak 
anomaly, and only at a single spectral point (the band center). 
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